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ABSTRACT 


Assume a technology that permits undistorted laser beam 
propagation from the aft section of a streamlined turret. 
A comparison of power on a distant airborne target is made 
between a Single aperture in a large scale streamlined turret 
with a turbulent boundary layer and various arrays of aper- 
tures in small scale streamlined turrets with laminar flow. 
The array performance is mainly limited by the size of each 
aperture. From an array one might expect, at best, about 
40 percent as much power on the target as from a single 
aperture with equal area. Since the turbulent boundary layer 
@n the large single-turret has negligible effect on beam 
quality, the array would be preferred (if all development 
efforts were esentially equal) only if a laminar wake is an 


@perational requirement. 
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Pr oteOr SYMBOLS 


a energy containing turbulence wavelength, m 


2 


> 
> 


incremental area element in the target plane, m 


A @iemarealwithin “the Bucket", m- 

B area of an aperture, m’ 

oF Cometalte ome roportionality, watts m 

Se constant of proportionality, m™ 

C. lift coefficient 

d interaperture spacing, m 

D. Single turret aperture diameter, m 

Da individual aperture diameter in the array, m 

F Dita tac cor 

h altitude, km 

H scale height, km 

I irddvanee on the target, watts m ~~“ 

Da irradiance on axis on the target for a single aperture, 
watts m-? 

1a aperture factor 

Ie, Eadtianmtmemittance of the aperture, watts m * 

i; irradiance at a point in the target plane at grid 
coordinates 1, j, watts m 

J, Bessel function of first order 

"4 wavenumber, m 

I aerodynamic lift force, 1lbf; boundary layer thickness, m 





v9) 


a = WM Fors 
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me) 


maximum array dimensions, m 


the number of grid points within "the Bucket"; Mach 
number 


the number of apertures in an array 
power projected, watts 


perpendicular distance between aperture and target 
planes, m 


Reynolds number 

wing surface area, ft? 
Pence nmareline. Kk 

Stagnation temperature, K 
flight velocity, m sec.’ and knots 


distance from leading edge of flat plate, m 


extinction coefficient for coherent beam in the 
turbulent boundary layer, m™ 


ratio of heat capacities 

geometric path difference between a distance extending 
from a point in the target plane to a reference 
position in the j® aperture, and the distance R, m 


angular divergence from beam axis 


proportionality constant which relates refractive 
index to air density 


radiation wavelength, m 
turbulence integral scale, m 
kinematic viscosity, m> sec 
dynamic viscosity, kg msec” 


atmospheric density, kg m~° 


local change in density, kg m~ 
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I. INTRODUCTION 


Propagation of laser energy from an aircraft turret 
designed with today's technology to an airborne target is 
unacceptably degraded when the target is in the quadrant 
between 1209 - 240° relative bearing due to the turbulent 
boundary layer and separated flow on the aft section of the 
turret. For Surface to Air Missile (SAM) defense, this is 
an especially interesting problem since the most dangerous 
threat to an aircraft is likely to approach from behind. 
meee la is a diagram of such a turret. The fairing is 
needed to curtail heavy oscillatory aerodynamic loading. 

Were the flow around such a turret laminar, an undis- 
meee Deam could be projected into the volume which the 
Separated flow presently inhibits. Reduction of turret scale 
could - if sufficient - result in laminar flow. However, at 
Meeecatt flight velocities and turret scales of practical 
size, a separated flow region will exist in any case 
eeue. lb). 

A future development may permit undistorted beam propa- 
gation from a window on the aft section of a streamlined 
body. In this event, a similar, though less severe beam 
degradation problem would probably still exist on the large 


scale due to the turbulent boundary layer (Fig. lc). Note 





the large scale turret would have a turbulent wake. A laser 
beam propagated through the turbulent wake would be severly 
degraded. 

The small scale streamlined body however, could achieve 
laminar flow (Fig. ld). Such a body with a future technology 
window could project a beam without degradation from separated 
regions or turbulent boundary layers. Beam dispersion in a 
laminar wake is not significant. Reduction of body size 
would however require more apertures to project an amount of 
power equivalent to that from a single larger aperture. An 
array of many such small turrets in laminar flows with total 
aperture area equal to the area of the single aperture may 
be distributed around the aircraft away from aircraft turbulent 
wakes. If properly configured, such an array will achieve 
greater power on a target aft of the aircraft than any 
other of today's technology turret designs, and perhaps 
greater power on a ee than the tomorrow's technology 
Single-turret design. Henceforth, discussion is limited to 


laser turrets designed using this future "window" technology. 





TODAY'S TECHNOLOGY (blunt base bodies on fuselage) 
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ice ee eoeN! AND FUTURE TECHNOLOGY TURRET DESIGNS 
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eee EOD 


One seeks a quantitative comparison between the per- 
formance of the single larger aperture projecting a beam 
through a turbulent medium and an array of small apertures 
projecting beams through laminar flows which converge on a 
Mistant airborne target. For identification, the Single 
Aperture projecting through a Turbulent boundary layer will 
be SAT; the Array of apertures projecting through Laminar 
flows will be AL. 

Isolation of the effects of the turbulent boundary layer 
on the single aperture and of the performance of the array 
in laminar flows is required to determine if greater power 
on the target can be achieved by the array. To do so requires 
the following assumptions: 

1. No power losses in the atmosphere. 

2. No thermal blooming effects. 

eNO )I1tCter. 

4. A perfect pointer/tracker. 

A reference area is defined as that area contained by 
the first zero of irradiance of the Airy disk which would be 
produced on the target in vacuo by the single aperture. 

This reference area is chosen believing that a single aper- 


ture design delivers sufficient destructive energy to a 


JU) 





target when projecting through a laminar flow and that energy 
outside of this area will not contribute to target destruc- 
tion or will be insufficient to disable it. This reference 
area will be called "the Bucket". 

Powerman the Bucket (PIB) tor the SAT and the AL will be 
compared. If the PIB for the array is found to be greater 
than that for the SAT, then the array would perhaps be more 


effective against targets approaching an aircraft from behind. 


ide) 





Pee iemienc LON OF LASER BEAMS 


The distribution of monochromatic (laser) irradiance 


Sia target is - among other things - the result of diffraction. 


Ewer SINGLE APERTURE 

For the single circular aperture in the Fraunhofer region 
without any dispersive atmospheric effects the resultant 
distribution of irradiance is the well known Airy disk. 
Calculations show that 83.8% of the incident power arrives 
within the area contained by the first zero of the Airy 


e@mstribution [Ref. 1]. 


Mee ARRAY OF APERTURES 

The array possesses maximum dimensions which are larger 
than the single aperture. Thus, the Fraunhofer condition 
may not be satisfied for the array at a target distance 
where it is satisfied for the single large aperture. 

Figure 2 shows an array of apertures where 0, is the 
source plane, 0 is the aperture plane, and 0, is the target 


plane. 


Ny 








Fig. 2: COORDINATE SYSTEM FOR AN ARRAY DIFRACTION PROBLEM. 
REEROBUCED FROM STONE. LREF. 2] 


Fomeaeplane wave incident on the aperture plane the 
Fraunhofer condition is satisfied when the maximum dimension 
Of the array creates a phase difference less than 0.1A in 
the target plane. This occurs when either of the following 
Semaitions are met. 

1. The beam is collimated, and L,/RA 221.4 

2. The array is focused at the target plane. 

A lens has the effect of bringing the Fraunhofer region 
closer to the aperture plane. It does this by lengthening 
miemceometric path length of light passing through the 
Memeral part of the lens such that the Fraunhofer condition 
1s satisfied near the focal point of the lens. Phase retard- 
ation across the aperture plane will have the same effect. 
This retardation can hopefully be accomplished by adaptive 
optics so that targets at close distances will be in the 


maaunhoter region. 
18 





For simulation. on a computer however, it is-much-easier 
to chose R>>Le/a than. to’: insert phase retardation across the 


apertures. 


C. IRRADIANCE ON A TARGET FROM AN ARRAY OF APERTURES 


The irradiance in the target plane is 


I.F (1) 


where Ih is the aperture factor, F is the array factor, and 


C, is a constant of proportionality. I, is a function of 


1 


aperture size and shape, distance to and location on the 


target, and the wavelength. For identical circular apertures 


2 
I = Agen or) (2) 
where 
3 (3) 


Feethe array factor, is the product (Ref. 3] of the 
conjugates of the sum of the complex phases at each point. 
The phases involved are a function of the array geometry, 
distance to and location on the target, and the wavelength. 


N ica? 
cS 


"TI 
ui 


(4) 
) 


where 


a = -xsin 9, San >, iS) 


1 





Refer to Figure 2 for the geometry and definition of 0, and 
>: C is chosen so that the irradiance on axis in the 
target plane is N* and has the units of watts/m*. Substi- 
tuting (2) and (4) into (1), the irradiance becomes 
Ps fon) ixa 6)” 
I= C a xen) || 5 ye (6) 

To visualize the distribution of irradiance at the target 
and to calculate PIB, a grid of points (X,5° Y94) is created 
in the target plane (see Fig. 2). At each point the phases 
and the argument of J, are computed. This yields vr at 
mac point; li is stored in an array. Those grid points 
within "the Bucket" contribute to the summation of PIB in 


the following manner: 


The PIB for a continuous intensity distribution would be 


ay I(x,y) dxdy (7) 


bucket 
This 1S approximated numerically as 


bucket radius 
PUB wy I,, AA (3) 
Ly 


which becomes 


bucket radius (9) 


A 
PD aay 2 I 
1 3 
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if the following further approximation is made: 
JAA == (10) 


mas approximation of AA improves as the density of points 
in the grid becomes greater. Calculation of PIB for the 
Airy disk by this method with about 80,000 points in the 
Pmemet yields 85.5% of the total Power Projected (P in 
the bucket. The published value is 83.8%. The program is 
in Appendix A. 

The PIB obtained from Equation (9) 1s compared to the 
meal power projected by the array to find the Fractional 
fetemern the Bucket (FPTB). FPIB is the true measure of 
the "goodness" of a system of one or more apertures. 

Eis 


AL 
FPIB, = + Gu 


Pp 
Power projected is found by use of the relation for the on 





axis intensity from a single circular aperture. 


2 
i Be (12) 
O 2 AR 


eememio is the intensity of the field in the aperture, B 


is the aperture area, A is the wavelength, and R is the 


distance to the target from the aperture(s). So that for 
an array 
2 
Z NBs , 
to 7 4, XR (13) 


AA 





The total power projected by the array is 


Ss. we (14) 


Radio astronomers use arrays of antennas; several papers 
[Refs. 4, 5] on radio astronomy were examined for concepts. 
Further, insects have eyes which contain a large number of 
elements in an array. The principals of insect vision [Refs. 
6, 7] were also investigated. 

The irradiance on the target plane was calculated for 
several array configurations; conditions for Fraunhofer 
diffraction existed. Figure 3 shows the FPIB achieved by 
each configuration. The best FPIB achieved for a 10-cm 
aperture was 48%. 

For arrays in the triangular configuration (probable 
candidate for use on aircraft) with array aperture area equal 
to the area of the single aperture, the ratio D,/D. will 
also determine FPIB. Figure 5 shows this trend which is very 
insensitive to interaperture spacing so long as the Fraunhofer 
condition is satisfied; the behavior is true in general for 
eum@earray Configurations and is qualitatively apparent from 


Paoure 4. 
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EPIBS 35% FPIB=34% 
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Fig. 3: FPIB FOR VARIOUS ARRAY CONFIGURATIONS 
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CIRCULAR 9-ELEMENT ARBAY TRIANGULAR 3-ELEMENT ARRAY 


m d=1.5D, 


FPIB=482% FPIB=68% 


TRIANGULAR 9-ELEMENT ARRAY TRIANGULAR 15-ELEMENT ARRAY 


q2 0. BR wor. 07 7 5 m d=0.2 m 
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meee 4a: LRRADTANCE ON THE TARGET FROM A SINGLE 
CEREUEAR APERIURE 
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Fig. 4b: IRRADIANCE ON THE TARGET FROM A TRIANGULAR 
5-ELEMENT ARRAY 
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miceeere: LRRADTANCE ON THE TARGET FROM A CIRCULAR 
9-ELEMENT ARRAY 
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Fig. 4d: IRRADIANCE ON THE TARGET FROM A TRIANGULAR 
9-ELEMENT ARRAY 
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ieee DY NAMECS RELATED TO TURRETS 


A. SEPARATED FLOW REGIONS 

Behind blunt bodies even at moderate Reynolds numbers 
turbulent separated flow will occur. Thus, the single large 
aperture in a circular turret must project through a region 


of turbulent separated flow such as that shown in Figure la. 


B. WAKES FROM TURBULENT BOUNDARY LAYERS 

Regardless of the body geometry, when the Reynolds number 
is greater than the critical value, a turbulent wake will be 
present behind the body. Propagation through a turbulent 
wake degrades beam quality more than either a separated flow 


or a turbulent boundary layer. 


fee TRANSITION TO TURBULENCE IN THE BOUNDARY LAYER 
An array of small scale streamlined turrets can maintain 

laminar flow over the many smaller turrets if the geometric 
scale is small enough. An upper bound on the irradiance 
which mirrorscan withstand before distortion and failure 
requires an equal area in the array as the single aperture 
to project equal power. The number of array elements then 
depends on the maximum aperture size that still maintains 
laminar flow since it is determined that arrays of large 


apertures achieve greater FPIB than arrays of small apertures. 





The array apertures are located inside of aerodynamically 
streamlined bodies with X,/D,=6 * (see Fig. ld). This design 
is proposed because of transition at a high Reynolds number, 


somewhere above R=2x10° Ree. 8]. 


(a = RY Ca) 
ie 6 D, 


Use of equation (15) and typical values of dynamic vis- 
cosity in the atmosphere as a function of altitude [Table l, 
Ref. 9] yields a critical velocity below which there is 


laminar flow over the turret in the array. 


Wap iE 1 
MevieEMat iG VisCOSITY VS, ALTITUDE IN THE ATMOSPHERE 


height (km) vx10) _ {m'/sec) 
. 2 
.5841 
an 
. 8648 
.0289 
r2ato 
.4164 
6456 
.9032 
.1936 
een 
.8948 
.5601 
3389 
.2508 
.3185 
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*For simplicity in these theoretical analyses it is assumed 
that the streamlined turret diameter is equal to the aperture 


(mirror) diameter in Figures lc syeval lial 





For various aperture geometric scales, Figure 6 shows 


@ritical velocity at altitudes from sea level to 15 km. 









ALTITUDE (KM) 
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[he /- ee 
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MERCRART VELOCITY (KNOTS) 


Peace Orme nt leat VELOCITY VS. ALTITUDE 


Figures 5 and 6 reveal that aperture geometric scale, the 
altitude where laminar flow is maintained, and FPIB for the 
array are all closely linked. It should be noted that an 
aircraft capable of transporting the laser system may not 
be able to fly at the velocity which corresponds to certain 
altitudes and aperture geometric scales. For example, if 
55% FPIB is desired from a triangular array, this requires 


10-cm apertures in the array if D_=50 cm. However, at an 


By 





Peestnae Or 5 Km transition from laminar to turbulent flow 
occurs at about 145 knots for a 10-cm turret. A large-body 
Hote dircratt probably cannot fly at this speed at an altitude 
@z 5 km. 

Thus, an estimate of regions of aircraft operability is 


needed; recalling the expression for aerodynamic lift: 


L = 'spV*SC, (16) 
Then 
ZG 3 
v= set 7 fale 


Zmemetnie Minimum velocity for an aircraft in a no-stall condi- 
mwome FOr most large-body jet aircraft L/S is typically 
about 100 lbs per square foot, and a practical maximum for 
eeepc assumed to be 1.2. Ve as a function of altitude has 


L 
been computed and is also shown in Figure 6. 
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enue GE OF SIURBULENT BOUNDARY LAYER 
ON POWER IN THE BUCKET 
A. SELECTION OF MODEL 
An estimate of PIB for the large single aperture is sought. 
The literature contains theoretical discussion of beam quality 
for propagation through turbulent boundary layers of uniform 
thickness [Ref. 10]. The angular dependence of beam intensity 


for the SAT is estimated by such a model. 


feeecoliIMATION OF PIB FOR SAT 

Sutton has calculated the time averaged angular dependence 
of beam irradiance for several cases of beam diameters (D.); 
energy containing turbulence wavelength (a) and boundary 
layer thickness (L) for the above model. 

Data points were taken from Figure 2 of Sutton [Ref. 11] 
Seems plots T/i,vs. 9D./A. The trapezoidal rule was used to 


meMmerically integrate 


1622 
a Gis) 
PIB. ay | T. Oe BXal 6 (sy) 
| oO 
where 
_ 98Ds (19) 


r 





A number representing PIB is obtained for each case of Di/a 
and aL.“ The PIB obtained in this manner is converted to 
FPIB by comparison to the result of the above integration 
using IE (EJ 7 for the Airy distribution, which is known. 


PIB 
SAT 
Pope = eg 3. 84 (20) 
ooh 7 ENO eae 


Table 2 shows these results. The program is in Appendix B. 


TABLE 2 
FPIB FOR SELECTED VALUES OF D./a AND aL 


aL 


100 





“a is nearly equal to A, the turbulence integral scale 


BRef. 12] 


Al 
Gt 





Lm indecheomnceoiem ao: fable 2 which most closely corres- 
ponds to the case of the single aperture streamlined turret, 
the following equations are used to estimate aL and Di/a. 

The equation for the extinction coefficient in the turbulent 


boundary layer is obtained from Sutton [Ref. 15]. 
eek oan oN CZ) 


where k is the wavenumber, An the change in refractive index, 
and A the turbulence integral scale which Sutton [Ref. 14] 


estimates as 


A = 0.1L (22) 


L is the path length through the turbulent boundary layer or 
the boundary layer thickness. This length is estimated using 
the equation for the flow of an incompressible fluid in a 
turbulent boundary layer [Ref. 15] beginning at the front 


edge of a flat plate. 


: 6.2 
fees 0.37 ux (235) 
oV 


The variables u and o have the usual meanings, V is the air- 
craft speed, and X is the distance from the leading edge 
(front of the streamlined body in this case}. The beam 
emerges from the rear of the body at a distance perhaps 4/5 
famene body length back from the front, or 5D. <A reasonable 


estimate for L is therefore 





02 
2 0+8 
L = 0.37 ba D. (24) 


Let C, = 0.37. (625u/pV)%*, which is clearly dependent on 


auGatmd@de and velocity. Substitution shows that 


mec (hen) DL" (25) 


It is known that 


an = K(®o) 52 (26) 
Or 
-h/H Ae 
An = Koy e — (ao) 


where K is the constant in the equation relating refractive 
index to density, i.e., n=l + KP/9,, H is the scale height, 
Ap is the change in local density, and p, is the ambient elec 


density. From Sutton [Ref. 16] one obtains themreomation f£or the 


Mieeeterm in Equation (27). 


a eel ( ae ou | (28) 


om = 


O 


where Oe is the density of air on the surface of the turret. 


Assuming an ideal gas and constant pressure across the boun- 


dary layer 


Ap . go. (1 - (29) 
Pe ee 


Gl 
~) 





Liepmann and Roshko [Ref. 17] provide an expression for the 
stagnation temperature ue = ieee as a function of the ambient 
temperature (Two), M the Mach number, and y the ratio of 


heat capacities. 
is =e ( A 1+ «| (305) 


Samce Stagnation temperature is equated to wall temperature, 
Pamation (50) is valid only for Prandtl number of unity. 
From the above equations it can now be shown that 

ier ns 2 


aL = 2@%)*10-°K, e7 2 h/H 3 BY is 


ee Y=] M2 2. S (31) 
i 


iveical values for substitution into (31) are H = 7.6 km, 
Gees km, D.= 0.3m, y = 1.4, K,, = 5Sxl107', M = 0.3, and 
c (5 km, 250 knts) = 0.16 m°*. These values produce an 
extinction number (aL) of about 1.1x10 °*. From Figure 2 of 
Sutton [Ref. 18] it is evident that the beam is negligibly 
affected by the turbulent boundary layer and therefore 
very close to the diffraction-limited case regardless of 


the ratio Di/a. 





VI. CONCLUSIONS 


Laser turrets designed with today's technology may be 
ineffective against an airborne target aft of the aircraft. 
An array of apertures using a future technology window and 
bmoper phase retardation can achieve reasonable Fractional 
mewer in the Bucket (FPIB) on a target approaching an aircraft 
meom behind. This FPIB is largely dependent on the aperture 
diameter of the array elements and is relatively insensitive 
to interaperture spacing as long as the Fraunhofer condition 
is satisfied. 

However, possession of the technology for a small scale 
window (see Fig. ld) for use in the array would probably also 
yield the technology for the large window (see Fig. lc) soon 
thereafter. Given the flexibility of both large and small 
window technologies, the single turret shown in Figure lc 
is by far the better. 

Even though the array achieves laminar flow and small 
windows may be easier to develop, such a system incurs the 
Meee Of a collosal fire control problem, phase control 
problem, and the operational restriction of flight velocities 
near the aircraft stall speed. Surprisingly, the single 
turret design is negligibly affected by the turbulent boundary 
layer and therefore nearly diffraction limited upon exit from 


the boundary layer. Though it does create a wake which is 


aS 





impwrene, Geasomaple advance detection of targets approaching 
from behind should avoid the necessity to try to project 


directly aft through the wake. 
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APPENDIX A 
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APPENDIX B 
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TRAPEZOIDAL RULE 
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APPENDIX B (continued) 
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